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Abstract 

Finite temperature Euclidean SU(2) lattice gauge fields close to the deconfine- 
ment phase transition are subjected to cooling. We find relatively stable or abso- 
lutely stable configurations with an action below the one-instanton action Si ns t = 
2tt 2 both in the deconfinement and the confinement phases. In this paper we at- 
tempt an interpretation of these lowest action configurations. Their action is purely 
magnetic and amounts to S/Si ns t ~ N t /N s , where N t (N s ) is the timelike (space- 
like) lattice size, while the topological charge vanishes. In the confined phase part 
of the corresponding lattice configurations turns out to be absolutely stable with 
respect to the cooling process in which case Abelian projection reveals a homoge- 
neous, purely Abelian magnetic field closed over the "boundary" in one of the spatial 
directions. Referring to the dyonic structure established for the confinement phase 
near T c and based on the observation made for this phase that such events below 
the instanton action St ns t emerge from dyon-antidyon annihilation, the question of 
stability (metastability) is discussed for both phases. The hypothetically different 
dyonic structure of the deconfinement phase, inaccessible by cooling, could explain 
the metastability. 



1 Introduction 

In the confined phase, below but sufficiently close to T c , applying cooling to Monte Carlo 
generated SU(2) lattice gauge fields has shown a dyonic structure of metastable action 
plateaus [H El El Ej ■ The dyons themselves correspond, in a good approximation, to the 
constituents of Kraan-van Baal-Lee-Lu (KvBLL) caloron solutions jSl El U\- The same 
cooling applied to lattice configurations in the deconfined phase acts entirely differently. 
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There is no remnant dyonic structure. Instead of this, there are metastable events on some 
lowest action plateau (actually significantly below the one-instanton action S inst = 2n 2 ). 
These configurations have vanishing topological density, and the action is purely magnetic. 
Such configurations have been already observed and discussed many years ago in papers 
by Laursen and Schierholz jS], and Veselov and Polikarpov 0. 

In the present paper, from the perspective gained with the dyonic structure at high 
enough temperature, we try to give a new interpretation of the lowest action configura- 
tions seen in Refs. [HI E] • For this goal we study them both in confined and deconfined 
phases. In the confinement phase we are in the fortunate position that we can observe 
the parent configurations which sometimes evolve into the configurations considered here. 
In all observed cases this parent configuration was a dyon-antidyon pair. Contrary to 
this, in the deconfined phase the cooling technique has been unable to exhibit potential 
parent configurations in the form of action plateaus, i.e. approximate solutions of the 
lattice equations of motion (metastable plateaus of action). We have critically examined 
the 't Hooft-Polyakov monopole structure of the cooled configurations emerging from the 
deconfinement phase that was suggested in Ref. [H]. We do confirm to see a minor fraction 
of magnetic configurations which resemble 't Hooft-Polyakov monopoles at a first view. 
However, even for these rare events we find only a weak correlation between the local- 
ization of magnetic action and the positions of monopoles defined either in terms of the 
Polyakov line or of the magnetic charge in the Weyl gauge. 

On the contrary, the similarity of the action dependence for these configurations on 
the spatial size of the lattice suggestively points toward their common nature. Most likely, 
in both cases we would expect quantized magnetic fluxes. The returned flux (unavoidable 
in the maximal Abelian projection for periodic boundary conditions and manifesting itself 
as a Dirac string) can be visualized as Dirac sheet (swept out by a Dirac string moving in 
Euclidean time). Hence, we adopt the name Dirac Sheets (DS) for this class of observed 
cooled lattice configurations. DS configurations are known as exact solutions of the lattice 
field equations in U{\) LGT [TO] . 

The paper is organized as follows. 

In Section |21 we will provide all necessary lattice definitions, in particular the observ- 
ables considered in order to identify KvB and DS solutions. 

In Section |21 we report on the statistics and the properties of DS events observed both 
in the confined and the deconfined phases. 

Section |U contains our conclusions. 

2 Production and characterization of DS solutions 

Throughout this paper SU{2) gauge theory in four-dimensional Euclidean space is con- 
sidered on an asymmetric lattice with periodic boundary conditions in all four directions. 
The respective ensembles of configurations have been created by heat bath Monte Carlo 



2 



using the standard Wilson plaquette action 



(1) 





with inverse coupling f3 = For simplicity the lattice spacing is set equal to a = 1. 

The lattice size was N% x N t with the spatial extension N s = 8, 10, 12, 16, 20 and with the 
inverse physical temperature T _1 = N t — 4. For N t = 4 the model is known to undergo 
the deconfinement phase transition at the critical coupling f3 c ~ 2.299 [H]. Throughout 
this paper we will use two ensembles with (3 = 2.2 < (3 C and /3 = 2.4 > (3 C . 

The equilibrium field configurations in both ensembles have been cooled by iterative 
minimization of the Wilson action S. In one or another form, cooling is used in order to 
smooth out short-range fluctuations, while (initially) leaving untouched some large-scale 
properties of the configurations. The cooling method applied here is the standard relax- 
ation method described long time ago in [12] and was used for investigation of instantons 



This method, if applied without any further limitation, easily finds approximate so- 
lutions of the lattice field equations as shoulders (plateaus) of the action as a function of 
cooling steps (relaxation history). Here we shall concentrate on smoothed fields at the 
very last stages of cooling, using a stopping criterium which selects the plateaus in the 
interval of action S < 0.6 S inst .. 

The emerging gauge field configurations were analyzed according to the spatial distri- 
butions of the following observables: 

• action density computed from the local plaquette values: 



(see eq.flU)); 

• topological density computed with the standard twisted plaquette discretization: 




(2) 




(3) 



• Polyakov loop defined as: 




(4) 



with 




(5) 



t=i 



where the Ug^ represent the links in time direction; 
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non-stationarity defined as: 

6 t = y %2\s(£,t + l)-s(£,t)\/S; (6) 



x.t 



violation of equations of motion A 

A = i^r E l tT IV** - u»*)(Uw - u x ,,Y] , (7) 



where 



U x ,n — ^U XjV U x+f>)IJ ,Ul + ~ jV + Ul_~^U x ^o^U x+ ^y^ (8) 

is the local link x, \i being the solution of the lattice equation of motion, with all 
degrees of freedom coupled to it held fixed. The factor c is just a normalization of 
the staple sum such that U Xtll G SU{2) x . 

Abelian magnetic fluxes and monopole charges defined within the Weyl (or gen- 
eralized temporal) gauge, doA = 0, and the maximally Abelian gauge (MAG). 
The first one is obtained from the Polyakov gauge (PolG, achieved by diagonalizing 
P(x) in each lattice site, followed by Abelian gauge transformations which render 
the (then diagonal) temporal links U^ T -,4 independent on time r. The latter is found 
by maximizing the gauge functional A 

A[g] = l E tr (^V3^V3)> ( 9 ) 



under gauge transformations U x ^ — > C/|„ = g{x)U x ^g\x + fi). The abelianicity is 
the maximum value of this quantity divided by the number of links. 

In both cases, Abelian link angles 9 X)ft are then defined by Abelian projection onto 
the diagonal U(l) part of the link variables U Xjfi € SU(2). According to the 
DeGrand-Toussaint prescription a gauge invariant magnetic flux B p through 
an oriented plaquette p = (x,/iz/) is defined by splitting the plaquette 6 P = 9 x>fi + 
x +fi,is-0x+u4i-0x,v into <d p = Q p + 27m p , n p = 0, ±1, ±2 such that Q p G (-7T, +7r]. 
The magnetic charge of an elementary 3-cube c is then m c = ^ ^2 p& q c @p • 

For the cooling procedure of equilibrium gauge field configurations we have kept the 
standard periodic boundary conditions on the AD torus. 

Finally, cooling was stopped at some (n-th) cooling iteration step when the following 
criteria for the action S n were simultaneously fulfilled: 



S„ < 0.6 S: 



inst j 



lr rhe replacement U Xtfi —> U XilJl is exactly the local cooling step as applied throughout this paper. 
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• S n — 2 o n _i + S'„_2 < 0. 

The last condition means that the relaxation just passed a point of inflection. As we 
have empirically observed, the point of inflection always coincides, within an accuracy of 
plus/minus one global cooling step, with a minimal violation of the equations of motion 
A. This can be understood as follows. If violation of equations of motion A = 0, i.e. 
equations of motion are fulfilled we are in the local minimum of the action where its 
variation is zero. If violation of equations of motion A has the minimum the variation of 
the action has also the minimum and second variation of the action is zero what means 
that the action goes through the point of inflection. 

3 Properties of Dirac Sheets in the confined and the 
deconfined phases 

3.1 Cooling 

We have investigated DS events on lattices iV s 3 x N t with N t = 4 and N s = 8, 10, 12, 16, 20. 
The statistics and the mean actions of various types of configurations selected by the 
cooling process are presented in Table 1 and Figures ^ and El where the dependence of 
S I Si ns t on N s /N t and N t /N s is shown to have the tendency Sos/Sinst — > N t /N s . Their 
properties are summarized in Table 2 and illustrated by the remaining Figures. We have 
found DS becoming very stable at action values ~ Si nst ■ N t /N s . The (color-)electric 
contribution to the total action is very small compared with the magnetic contribution. 
Moreover, they are perfectly static with the values of nonstaticities 8% shown in Table 
2. Employing MAG we have convinced ourselves that they are almost Abelian (see the 
abelianicity A in Table 2), and only a small fraction contains MAG monopoles. 

In the confinement phase it happens quite rarely that they appear directly in the result 
of the cooling process. In all cases observed they appear after dyon-antidyon pairs have 
been observed at S m S- mst which annihilate in the final stage of relaxation. The Abelian 
monopole content of DS in the deconfinement phase, if it is obtained by Abelian projection 
in the Weyl gauge, amounts to monopole-antimonopole pairs being present in 60 -j- 90% 
of those cooled configurations. 

We originally had found this type of solution for fixed holonomy boundary conditions 
(f.h.b.c.) 121 El El- These DS for f.h.b.c. were seen to be oriented very exactly in plane 
and to have non-zero action for plaquettes in one of the space-space coordinate planes 
(x,y), (x,z) or (y,z). They had the same action values S/Si nst as we have found lateron 
for the case of periodic boundary conditions. We had also found configurations which 
contained two DS orthogonal to each other at action values twice as large than for one 
DS. We shall not further comment on such events in this paper. 
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Table 1: The statistics of DS events on the lattices iV 3 x N t with N s = 8, 10, 12, 16, 20 
for (3 = 2.20 (confined phase, upper row) and j3 = 2.40 (deconfined phase, lower row). S 
denotes the average action of the observed events, SS the variance. 



p 


N./N t 


^ event 


S / Si ns t 


5S / Sinst 


2.2 


2.0 


33 


0.479 


0.007 


2.4 


2.0 


54 


0.471 


0.004 


2.2 


2.5 


18 


0.392 


0.009 


2.4 


2.5 


62 


0.387 


0.003 


2.2 


3.0 


12 


0.331 


0.002 


2.4 


3.0 


73 


0.323 


0.002 


2.2 


4.0 


11 


0.249 


0.0007 


2.4 


4.0 


66 


0.245 


0.002 


2.2 


5.0 


7 


0.200 


0.000 


2.4 


5.0 


28 


0.195 


0.002 



3.2 Can DS configurations be viewed as 't Hooft-Polyakov 
monopoles in the deconfinement phase ? 

In contrast to previous parlance (these configurations have been called "monopole" (M) 
configurations in Refs. [B1IH]) we have called them here Dirac Sheet (DS) configurations 
from the beginning. Now we want to provide some more facts supporting this inter- 
pretation. In the beginning let us recall several features which originally suggested the 
interpretation as 't Hooft-Polyakov monopoles. 
There should be 

• a minimum of the dynamically generated "Higgs field" tr(v4 /T) 2 

• a maximum of the (almost purely) magnetic action density trB 2 

• a pair of pointlike Abelian magnetic charges, one carrying action and the other being 
spurious. 

We have studied these signatures for our configurations cooled down below the one- 
instanton level starting from the deconfinement phase, i. e. from equilibrium Monte Carlo 
configurations on a 16 3 x 4 lattice at f3 — 2.4. All these configurations have one maximum 
of the magnetic action and a relatively shallow maximum of the Polyakov line. A subset 
of these configurations, considered in the Weyl or maximally Abelian gauge, respectively, 
has a pair of corresponding static magnetic charges. The relative frequencies actually to 
find the indicated Abelian monopole structure is given in Table 2. 



6 



Table 2: The properties of DS events on the lattices Ng x N t with N t = 4, N s = 
8, 10, 12, 16, 20. For each iV s the events are presented in three rows: unstable (upper 
row) and absolutely stable (middle row) DS in the confined phase for (3 = 2.20 and DS in 
deconfined phase for (3 = 2.40 (lower row). 



N s /N t 


frequency of 
DS events 


fraction of DS 
with monopoles 
in Weyl gauge 


fraction of DS 
with monopoles 
in MAG 


A 


A 


5 t 


2.0 
2.0 
2.0 


3.8% 
3.2% 
4.4% 


21.1% 
59.1% 


5.3% 
6.8% 


98.8% 
100% 
98.0% 


0.223E-05 
0.944E-09 
0.355E-05 


0.372E-03 
0.253E-06 
0.889E-03 


2.5 
2.5 
2.5 


3.0% 
4.8% 
11.3% 


41.7% 
79.4% 


8.3% 
2.9% 


98.9% 
100% 
98.8% 


0.122E-05 
0.167E-15 
0.107E-05 


0.426E-03 
0.245E-12 
0.512E-03 


3.0 
3.0 
3.0 


2.0% 
5.7% 
14.3% 


50.0% 
88.4% 




99.5% 
100% 
99.1% 


0.378E-06 
0.150E-15 
0.464E-06 


0.214E-03 
0.389E-12 
0.397E-03 


4.0 
4.0 
4.0 


0.3% 
5.0% 
16.0% 


72.9% 


14.6% 


99.9% 
100% 
99.5% 


0.385E-10 
0.760E-15 
0.665E-07 


0.103E-07 
0.108E-10 
0.149E-03 


5.0 
5.0 
5.0 


2.8% 
18.0% 


77.8% 


11.1% 


100% 
99.7% 


0.636E-16 
0.244E-07 


0.245E-11 
0.951E-04 



One of the cooled configurations with a Weyl gauge monopole pair is portrayed in 
Figure 01 where we show the profiles of action density and Polyakov line over the (x, z) 
plane which cuts the configuration and contains the maximum of action density. In Figure 
El the positions of the static monopole pair are projected onto that (x, z) plane. This event 
indeed shows a structure as reported in Ref. jHJ. 

For a more quantitative assessment, we have compared the positions of the maxima of 
three-dimensional action density and the Polyakov loop. We found the average distance 
equal to 7 ± 2 lattice spacings, whereas for the maximum of the action density and the 
position of the nearest Weyl gauge magnetic monopole the distances are somewhat less, 
5.2 ±1.6 lattice spacings. These distances are certainly smaller than the maximal three- 
dimensional distance 8a/3 ~ 13 but are difficult to reconcile with a static extended particle 
(monopole) interpretation. 

Monopoles obtained in the MAG are present only in approximately < 15% of all DS 
events (see also Table 2). One of the rare cases when a DS event in the deconfinement 
phase shows an Abelian monopole-antimonopole pair in MAG is presented in FigsEl and 
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3.3 More on the Dirac sheet interpretation 

At this point it might be useful to focus on the striking volume dependence of the action 
(common to both phases) and on the surprising stability of the cooled DS configurations 
(which distinguishes the confinement phase) in order to understand similarities and dif- 
ferences between the ways how these configurations originate from the respective vacuum 
of the two phases. 

The dependence of S/S inst on N s /N t and N t /N s is shown on Figures [T] and El In the 
confined phase for N s = 20 all 7 DS events were absolutely stable with Sps/Sinst = N t /N s . 
The number N t /N s can be understood as follows. Let there be a quantized homogeneous 
Abelian magnetic flux in some spatial (x, y or z) direction. The Abelian magnetic field 
is equal to = Ati/N]. Its action is equal to [TU] 

Sds = ~ 2 {l ~ cos(^))N*N t » ^- 2 {BlfNlN t = = S mst §- . (10) 

9o z z 9o 9o Jv s iV s 

All 7 DS events in the confined phase for N s = 20 (when put into MAG) show such 
an Abelian magnetic field. There arises the question: why for other (smaller) N s in 
the confined phase the fluxes are not always homogeneous and absolutely stable, and 
why the fluxes are unstable for all N s in the deconfined phase. The probable answer is 
that magnetic flux (during the process of DD annihilation) is not always closed over the 
"boundary" in some periodic spatial direction. If the size of dyons in the DD pair is small 
compared to the spatial size of the lattice, the annihilation is almost pointlike and the 
magnetic flux has a good chance to be closed. 

The size distribution of dyons depends on holonomy. In the confined phase the mea- 
sure of holonomy L = cos(27ro;) (oj being the holonomy parameter) is distributed in the 
neighbourhood of zero. For lower plateaus of action the distribution approaches more and 
more the semicircle law (Haar measure). Then uj is distributed over the range < uj < 1/2 
and the size of the lighter dyon as known from the KvB solution (Nt/Anu for < uj < 1/4 
and N t /4ir(l/2 — to) for 1/4 < uj < 1/2 measured in lattice spacings) varies from N t /ii to 
the maximal value possible on the finite lattice. In the deconfined phase, the holonomy 
becomes closer and closer to the trivial one (L w ±1) and the dyons (in the dyon-antidyon 
pair) are strongly delocalized. As mentioned in the Introduction, during cooling in the de- 
confined phase the dyon-antidyon pair itself does not become visible on a well-established 
plateau. 

The correlation between holonomy and stability of DS events is shown in Figure [7| for 
both phases. It can be understood if the configurations are really emerging from the an- 
nihilation of a dyon-antidyon pair. The Figure presents scatter plots where each DS event 
is represented by two points: (s min , holonomy) and (s max , holonomy) with s min and s max 
being the action density at sites where it is minimal and maximal, respectively). Pro- 
vided that the holonomy remains far enough from trivial, we obtain DS events from the 
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confined phase which consist of homogeneous Abelian magnetic fluxes. The homogene- 
ity is expressed by s m i n = s max and corresponds to the successful annihilation of more 
or less pointlike dyon-antidyon pairs. However, for values of holonomy close to trivial 
holonomy DS events in both confined and deconfined phase occur as unstable magnetic 
fluxes which are not closed as the result of annihilation of less localized (and "massless", 
i. e. low-action) dyons. The unstable DS in confined and deconfined phases have similar 
characteristics as can be seen from the first and third rows (shown separately for each 
N s ) in Table 2. 

So, unstable DS events in confined and deconfined phases are similar. There is no 
absolute gap between unstable and absolutely stable DS events in the confined phase. 
This can be an argument in favor of their common nature. The stable DS events found 
in the confined phase are purely Abelian magnetic fluxes. 

4 Conclusions 

We have generated SU (2) lattice gauge fields at non-zero temperature, both in the con- 
fined and the deconfined phases. We have cooled them in order to analyze the structure of 
the lowest action plateau (which in fact is below the one-instanton action). We have found 
certain structures ("Dirac sheets") that resemble homogeneous Abelian magnetic fluxes. 
The action dependence on the spatial lattice size N s favors such an interpretation. For 
the deconfinement phase, where an 't Hooft-Polyakov monopole interpretation has been 
advocated, the loose correlation between different possible definitions of how to localize 
the monopole as an extended heavy particle makes this picture less convincing. Instead, 
we have looked for another interpretation in terms of how the original dyonic structure 
(which is different in the two phases) becomes destroyed by the cooling process. In the 
infinite volume limit N s — > oo, these DS structures disappear. Therefore, we interprete 
them as artefacts of the finite lattice volume. 
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Figure 1: Action of DS events in the confined (triangles) and in the deconfined (circles) 
phases for lattices with N t = 4 and N s = 8, 10, 12, 16, 20 as function of N s /N t . 




Figure 2: Same as in Fig. Q]as function of N t /N< 
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Figure 3: Action density (left) and Polyakov line (right) of an DS event in the deconfined 
phase shown as function over the xy plane cutting through the maximum of action density. 



mon 




Figure 4: Abelian monopole charge density in Weyl gauge (summed over z) of the DS 
event presented on FigOl For static monopoles the monopole charge density is equal to 
±iVf = ±4. The smearing is due to the interpolation in MATHEMATICA graphics. 
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Figure 5: Action density (summed over t and z) of an DS event in the deconfined phase 
shown as function over the xy plane. 



mon 




Figure 6: Abelian monopole charge density in MAG (summed over z) of the DS event 
presented on Fig|S| For static monopoles the monopole charge density is equal to ±N t = 
±4. The smearing is due to the interpolation in MATHEMATICA graphics. 
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Figure 7: Correlation between the Polyakov loop L and the stability of DS configurations, 
illustrated by scatter plots showing, for each DS event, two points: (s min , L) (as triangles) 
and (s max , L) (as circles) with s min (s max ) being the minimal (maximal) action density 
of the DS configuration. The left coloumn presents DS events in the confined phase for 
lattices with N t — 4 and N s = 8, 12, 16, 20 (from up to down), the right coloumn presents 
DS events in the deconfined phase for the same lattice sizes. 
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